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We construct a class of microscopic electron models with exact SO (5) symmetry between antiferromagnetic and d 
wave superconducting ground states. There is an exact one-to-one correspondence between both single-particle and 
collective excitations in both phases. SO(5) symmetry breaking terms can be introduced and classified according to 
irreducible representations of the exact SO (5) algebra. The resulting phase diagram and collective modes are identical 
to that of the S0(5) nonlinear a model. 



PACS numbers: 74.20. Mn, 74.25. Ha, 71.10.-w 

One of the most interesting features of the high Tc su- 
perconductivity is the close proximity and the interplay 
between the antiferromagnetic (AF) and the d wave su- 
perconducting (dSC) phases. Recently, a theoretical for- 
malism was introduced based on the concept of a SO {5) 
symmetry between these two phases, and the resulting 
field-theoretical model describes the cuprate phase dia- 
gram and collective modes in a unified framework [Q. 
It was argued that the microscopic Hubbard model sup- 
ports an approximate SO{5) symmetry 

In this note, we construct a class of microscopic elec- 
tron models with exact iS'0(5) symmetry. In this model, 
degeneracy between the AF and dSC phases can be 
demonstrated exactly, and both the fermionic single- 
particle and the bosonic collective modes can be mapped 
onto each other with a precise one-to-one correspon- 
dence. This model can be used as a starting point around 
which iS'0(5) symmetry breaking interactions can be in- 
troduced and classified according to irreducible tensors of 
the SO{5) algebra. It is shown that the resulting phase 
diagram and the collective modes are similar to those ob- 
tained from the effective 5*0(5) nonlinear a model with 
anisotropic couplings . The purpose of this work is to 
demonstrate that the general 50(5) idea can be realized 
exactly by explicit microscopic Hamiltonians. The mi- 
croscopic information extracted from this class of mod- 
els, especially the behavior of the fermionic excitations 
across the AF/dSC transition would greatly complements 
the effective field theory approach. Within this class of 
models, we have a consistent microscopic theory of the 
AF/dSC transition. Since both the AF and the dSC 
states are stable infrared fixed points, it is plausible that 
one can deform the parameters so that the microscopic 
5*0(5) models can also serve as a paradigm for a much 
more general class of AF/dSC transitions, including those 
occuring in the high Tc cuprates and 2D organics. 

The first independent attempt to construct micro- 
scopic 5*0(5) invariant models was undertaken by 
Christopher Henley ^j. He independently made a cru- 
cial observation that if one replaces the standard d wave 
factor cospx — cospy by sgn(cosp3: — cospy), the 50(5) 
algebra introduced in ||^ closes exactly. 



It is easy to write down many electron models with ex- 
act SU (2) spin rotation invariance, because the electron 
operator Cpa forms a natural spinor representation of the 
SU{2) algebra. In writing down SU{2) invariant models, 
all we have to do is to contract the spinor indices in a 
natural way. Therefore, the first step towards construct- 
ing a 5*0(5) invariant electron model is to find a natural 
definition of a 50(5) spinor. Spinor representations of 
the 50(5) Lie algebra can be easily constructed using 
the Clifford algebra of five 4x4 Dirac matrices |^ satis- 
fying {F",r''} = 25ab {a,b= 1, ...,5), and the ten 5*0(5) 
rotation generators are given by F"'' = — i[r°, F**]. In this 
work we shall use the following explicit representation for 
the Clifford algebra: 
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where a = {ax, (Jy,<Jz) are the usual 2x2 Pauli matrices, 
and * means transposition. We define a 4 component 
spinor by 

= {cpT>Cpi,(/)x(p)cLp+Q p0^(p)cLp^g^^| (1) 

where (p) = sgn{cos Px— cos Py) = ±1, and Q = (7r,7r). 
Since we have two spin degrees of freedom at a given mo- 
mentum p, such a description must be redundant. In- 
deed, one can easily see that the spinors with momenta 
outside of the magnetic Brillouin zone is related to the 
spinors inside the magnetic Brillouin zone by an "R con- 
jugation" 



*P+Q = <?!)x(p)i?*! 



(2) 



the form R 



The R matrix is an invariant tensor of the 50(5) 
algebra enjoying the following properties: RV^R = 
-*F°, RV^R = *F°^ In our representation it takes 

or 

-1 

trix is related to the fact that the spinor representation 
of the 5*0(5) Lie algebra is pseudo-real. The dy matrix 
plays a similar role for 50(3)). The ^'p^ spinors obey 
the anticommutation relations: 



(The existence of such a ma- 
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(3) 



{*pa' *p'/3} = {*pa: *p'/3} = -'^-(P)^a/3 W',Q- (4) 

If we restrict p and p' to be both inside the mag- 
netic Brillouin zone, the right hand side of the second 
equation vanishes and the ^'p^ spinors commute in the 
same way as the Cp„ spinors. They can be used to 
construct the 50(5) vector order parameter and the 
symmetry generators: = |: I]p ^p^p^^^p+q, and 
Lab = ■|Ep*p^°''*P- Here Wp = w^p. Note that the 
definition of the tt operators (ii (2,3,4), ^-(2,3,4)5) differs 
from the ones used in previous works ||l|-|3|], where they 
are electron pair operators on the nearest neighbor (n.n.) 
sites. The problem with this kind of definition is that the 
commutator algebra does not close, and generates longer- 
ranged bonds. Naively, the condition for the closure of 
the 5*0(5) algebra appears to be over-constrained. The 
present work and Q start with electron pair operators 
with a long-ranged profile, given in real space by the lat- 
tice Fourier transform of 07r(p), 



(/)^(m,n) 
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(5) 



where R = (m, n) is a lattice point. It is truly remark- 
able that this simple choice closes the algebra exactly. 
Notice that while the tt operators have long-ranged pro- 
files, the dSC order parameter can still be short-ranged 
with suitable choices of Wp. Under the 50(5) rotations 
generated by the Lab, transforms as a proper 50(5) 
spinor 



[Lab,^pa]=-\iT''')a0^p0 



(6) 



for all values of p. Using these spinors, exact 50(5) 
invariant Hamiltonians can be constructed simply by 
proper contraction of the spinor indices. 

We start with the kinetic term, and write it as 



-^^kin — ^ ] ^p'^p,cr''P,' 
p,o- 



2 Z^^P^p^p- 
p 



(7) 



Ep, valid for n.n. 



We see that the property Ep+q = 
tight binding model, is crucial for this construction to 
work. In order to construct four-fermion interactions, we 
first note that a 50(5) spinor bilinear can in general be 
decomposed into a direct sum of a scalar, a vector, and an 
antisymmetric tensor, i.e. 4 x 4 = 1 + 5 + 10. Therefore, 
general 50(5) invariant four-fermion interactions can be 
expressed as 

J2 V, (p, p'; q)(vi/^r°*p+q)(*;,r'' v-q) 

p.p'.q 

+ E ^2(p,p';q)(*J,r'^''*p+q)(<,r'^''V-q) 

p,p',q 

+ E ^0(p,p';q)(^^*p+q)(*I,'V-q)- (8) 

p,p' q 



Since L„h(p,q) = ^'t^r'-^^p+q, na(p,q) = «'t,r»«'p+Q+q 
and p(p, q) = ^'j^^'p+q are the true 50(5) tensor, vec- 
tor, and scalar, respectively, for any p and q, their in- 
ner products naturally gives a manifestly 50(5) invariant 
Hamiltonian. 

Among three terms in iJjnt , we concentrate on the vec- 
tor interaction (first term) in all subsequent analysis, and 
assume a factorizable form Vi(p,p';q) = —Vi{c\}WpWpi] 
This form is not necessary, but simplifies calculations. In 
real space. 



i,n 



(9) 



Here, rfii and Ag are Neel and d-wave pairing order pa- 
rameters (operators) at site £ = Hg, but with extended 
internal structures determined by Wp. For the simplest 
choice Wp = 1, they become 
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)^(Rf - Rj)(c]^c]^ - cj^c]^). 
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Here we introduced two-component spinors tpi 

\c^,Cii) and xe = (- 

J2j fpTri'Ri — R-j 



) X \bif,bei) with bi^r = 
The pair wave function for dSC 
condensate is described by ip^ and is long-ranged. For 
the choice Wp = \ cospx — cos Py\, we obtain 



me 



E 0A/(R^ - - XpXe) (12) 
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where 



0Af(TO,n) 



1 + (-)™+'^ 



TT^ [(m -I- n)2 — 1] [(m — n)^ — l] 



(14) 



The interaction between centers of mass of m or A fields 
is controlled by Vi (R) . If we take Vi (q) to be a (5 function 
at q = Q, the A-part in iJint,i becomes the usual BCS 
reduced Hamiltonian for n.n. c?-wave pairing. If Vi (q) is 
taken to be a Lorcnzian around q = Q, the real space 
form of the spin interaction resembles the potential in- 
duced by the AF paramagnon exchange [^|-^ . 

It is not difficult to find degeneracy between AF and 
dSC states in the usual treatment of mean field theories. 
However, their excitation spectra are generally different, 
and quantum fluctuations may remove this degeneracy. 
In the 50(5) invariant models, symmetry not only en- 
sures exact degeneracy of the ground states, but also 
ensures exact one-to-one correspondence between their 
excitation spectra. This fact is formulated as follows: 
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Theorem 1: If Ivfo > is a ground state of a 5*0(5) 
invariant Hamiltonian with AF broken symmetry (say 
in the n2 direction), i.e. < \I'o|?T-a|*o >= S2,aA, then 
jv^Q >= e*^^"|^o > is a degenerate ground state 
with dSC broken symmetry (in ni direction), i.e. < 
'J'ol'^al^'o >= <5i,a^- Furthermore, aU excited states of 
the AF ground state can be mapped to excited states of 
the dSC ground state at the same energy by the e'^'^i^ 
operator. 

The proof of this theorem is elementary, since L12 com- 
mutes with the Hamihonian, and e~'^ ^^^nie^'S'^i^ _ j^^. 
In the following, we shall illustrate this powerful theorem 
in an explicit mean field calculation. We take a "gener- 
alized BCS reduced Hamiltonian" by selecting Vi(q) = 
^i<5q,Q in the vector interaction. The Green's function 
in the presence of a mean field (Up) = ■i(5'|jr°^'p+Q) is 
given by 

G^f,{p,p';cj) = -t I dte'^\T^p,^ml, f,{0)) (15) 
^ (c^ + £p)(5aff(5p,p> + A°r°^,5p,p.+Q 



Theorem 2: In a 50(5) invariant model, there are four 
gapless Goldstone modes associated with spontaneous 
symmetry breaking. 

In the AF phase, there are two spin wave modes and 
a TT doublet (4 = 2-1-2). In the dSC phase, there is one 
SC phase mode, and a tt triplet (4=1 + 3) 

As symmetry-breaking perturbations to the above 
50(5) invariant Hamiltonian, we consider two typical 
terms. One is the coupling to external fields Bat, 



a<b 



BahLn 



(16) 



A particular example of this field is the chemical potential 
_Bi5 = — 2/i, leading to = 2/iLi5. The other is the 
anisotropy energy 



E 

p,p',qa— 2,3,4 



J2 5(q)(*^f"*p+q)(<'r'^*p'-q) (17) 



where A" = -16ViWp 



between AF and dSC states. To study the spectrum of 
H = i/kin + -ffint.i + Hf, + Hg, we take g{q) = gSq^Q and 
again use mean field approximation. In the dSC phase, 
(rip) lies in the plane (7x1,715). We choose it in the ni 



This manifestly direction. Then the Green's function is given by 



50(5) invariant Green's function explicitly shows that 
the AF quasi-particles can be mapped onto dSC quasi- 
particles. In particular, the AF Green's function (in 
the 7i2 direction) can be obtained directly from the dSC 
Green's function (in the ni direction) by a simple rota- 
tion: G^^ -- 



(t^+ep-Ai)l(5p 



— i(T„ A„(5, 



P°p, p' + Q 



G^^(p,p',c.) = 



p"p,p' + Q 



(w + £p+Al)l(5 



If we take Wr 



1, 



the AF quasi-particles have a full s wave gap, while the 
dSC quasi-particles have a full d wave gap, with step 
discontinuity at (±7r/2, ±7r/2) points. For the choice of 
Wp = I cospx — cospyl, the dSC quasi-particles have the 
usual cospx — cos py gap behavior (Fig.l), while the AF 
quasi-particles have an anisotropic s wave gap with nodes 
at (±7r/2, ±7r/2) points (Fig. 2). Because the AF nodes 
are not "topological" , any small interactions will remove 
it ||ll|. In either case, the amplitude of the gaps are the 
same in both phases. 

In addition to the gapped single-particle excitations, 
the 50(5) invariant models also have four gapless Gold- 
stone modes. There are ten conserved Noether cur- 
rent density LJJ^(x) associated with the 50(5) symme- 
try, where /i (= 0,1,2) is the space-time index. In the 
small q limit, exact current conservation leads to the fol- 
lowing generalized Ward identity for the vertex function 
l^^'i.P + ItP) ip — P)) related to this current 

7;'(P + 9, P)) ^ = - J {^tGZl {P + l)- G-l {p)Tfp 

For a non-zero (n^), the rhs of the above equation has 
a finite 9 — > limit, therefore, the four vertex functions 
l^{p + QtP) with b ^ a must have gapless poles in this 
limit. This result is formulated in the following theorem: 



where Ap = -16VL(/'7r(p)wp ■'"k(nk) = Ao0,r(p)wp 
and the g term drops out because of symmetry mismatch. 

Ao is determined by the gap equation 1 = 16 Vi 

where i?k = \/ (ck — m)^ + ^k- ^^"^ choice Wp = 

I cospx — cospj,!, we have a usual cZ-wave gap. In the AF 
phase, {rip) lies in the (712,773,774) space. If we pick the 
714 direction, we have 



G 



AF 



(P,P',t^) = 

+Ap 



el-Al+i& 





"p.p'+Q 



(w_-|-£p)15 





,'+Ap 



-Al+iS 



Here w± = uzL^ and Ap = — 16 J^ki^^'^p'^^ + 5)(^*k) = 
TWpAo -|- Ag. Aq and A^ are determined by the gap equa- 
tion Ap = 16 X]k(^i^p'"'k + 9)^^- The real space form 
of Ap has an on-site contribution from A^ and a long- 
ranged contribution proportional to Ao0m(R-)- We see 
that the g terms leave dSC gap unaffected, while it re- 
moves the AF gap node (See Figs.l and 2). The ground 
state energy curves are shown in Fig. 3. The "superspin 
flop" transition from AF to dSC occurs at /i = for 
Hg = 0, while it occurs at a finite value of fic for 5 7^ 0. 
In this case, the AF/dSC transition is first order, with a 
finite jump in hole density Xc (See Figs. 3 and 4). 

While the above pictures are based on the mean 
field approximation, some exact statements can be made 
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about the AF/dSC transitions. 5*0(5) is a rank 2 al- 
gebra, and we can choose Q = — L15 and ^2 — — -Zj23 
as the members of the Cartan (maximal commutative) 
subalgebra. In addition, we have the Casimir operator 
C = X]a<6-^ab' which commutcs with all the generators 
and has eigenvalue l{l + 3). The set {Q,Sz,C) forms a 
Cartesian coordinate system labelling the quantum num- 
bers of all states in the Hilbert space. If we consider only 
states with even number of electrons, these states form a 
pyramid, with the I — singlet on top, the I — 1 vector 
next, and the I = 2 traceless symmetric tensor on the 3rd 
layer etc. States on the same layer are all connected by 
repeated actions of the 8 root generators. Therefore, 

Theorem 3: In SO{5) invariant models, it is sufficient 
to diagonalize the Hamiltonian at half-filling with Q = Q 
and Sz = 0. All the other states (with even number of 
electrons) in the Hilbert space can be obtained from these 
states through the action of {Sx, Sy, 7?, 7?^^). 

In this sense, states at half-filling fully determine the 
states away from half-filling. In the presence of the i?^ 
term, the vr^ and tTq, operators are exact eigenoperators 
of the Hamiltonian with eigenvalue ±2/i. Therefore, 7?^ 
commutes with the Casimir operator, and simply shifts 
the energy of the Q ^ states linearly without chang- 
ing the wave function of these states. In a system with 
spontaneous symmetry breaking, lowest states with dif- 
ferent / quantum numbers are separated by inverse sys- 
tem size. In the infinite system limit, these shifts due to 
chemical potential converge to the parabola as depicted 
in Fig. 3 

The symmetry-breaking terms, Hg and i/^ produce the 
mass gap in the Goldstone mode spectrum. For Hg = 0, 
the mass of the tt triplet mode is exactly 2|/i|. For finite 
Hg we employ the equations of motion (EOM) . 

As discussed earlier, Lats commute with i?kin + ^^int,i 
for arbitrary Vi(p,p';q). The only contributions to the 
EOM for Lab come from the 50(5) symmetry breaking 
terms, Hg and iJ^: 

[Lab, H] = —i ^(-Bac-^fac — BjjcLac) (18) 

c 

- i{Sa - St) ^ 5(q + Q){na(q), nf,(-q)}, 
q 

where 5a = 1 for a = 2, 3, 4, and = otherwise. This is 
precisely the same equation as that obtained in . EOM 
for rio is given by 

K,i7] = i^(ep+Q - ep)na{p,0) + i^Babnb (19) 

P b 

-iJ2iVM + Q) + ff(q + Q)6b) {nb(q), Lat(-q)}. 

where na(q) = | Ep"-a(p, q), ia6(q) = I Ep ^a6(p, q), 
and a form Vi(p,p';q) = —Vi{q) and the property 
^i(Q — q) = ^i(Q + q) has been assumed. 



Let us now consider the case where only the chemical 
potential Bi^ = — 2/x is introduced as external fields, and 
consider only the collective part in equation (|l^) ||l3|,|l^ . 
In the dSC state, we take (7ia(q)) ~ (ni)(5a.i<5q,o, and 
linearize eq.(p^ to obtain 

= 2l^i(ni)ii„ (a = 2, 3, 4) (20) 

n5-2/i = 2V^i(ni)ii5, (21) 

where Vi = Vi(Q). These equations should be compared 
with x'i'a — Lia with a = 2,3,4 and x("'5 ~ 2/x) = L15 
derived from the nonlinear cr- model ||l|. Equations (p^), 
(pT|), and (|2^) can be combined to give fis = and 
ha = 4(/^^ _ y^'^)na, where /^c = {ni)^/gVi and g = g[Q). 
Therefore, the energy of the triplet tt excitations in 
the dSC state is given by = 2v/ — /i^, which is 
also consistent with the result of [|l|. Similar calcula- 
tion in the AF state gives the energies of the tt doublet 
= 2(714)1/(7 (Vi + g)± 2//, where we assumed AF or- 
dering along 714. We therefore see that the two symmetry 
breaking terms g and /i partially compensate each other 
for the TT triplet and Q = —2 it doublet. 

In conclusion we have constructed a class of electron 
models with exact 5*0(5) symmetry. Both the fermionic 
single particle and bosonic collective modes of the AF 
and dSC phases are in one-to-one correspondence to each 
other. Energy levels are classified according to the 50(5) 
quantum numbers and the level crossing at the AF/dSC 
transition can be followed in detail. The fermionic single- 
particle gaps do not close at the AF/dSC transition, but 
instead rotate from one direction in superspin space to 
another. It is amusing to ask what experimental results 
this kind of ideal models would predict. The phase dia- 
gram would be identical to the one depicted in Fig. lA of 
reference with an insulating AF gap and a finite jump 
in chemical potential at half filling, phase separation or 
stripe ordering for doping range < 2; < a;c, a low energy 
spin triplet tt resonance in the dSC phase, and a "pseu- 
dogap" behavior in the high temperature phase above 
the bi-critical point, where the gap direction fluctuates 
(rotates) between AF and dSC character. 
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FIG. 3. The groundstate energy Eg in both AF and SC 
phases as functions of ^. 
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FIG. 4. The electron density x versus ^ in the presence of 
anisotropy energy Hg. 





FIG. 2. The antiferromagnetic (AF) gap. The solid (dot- 
ted) line is for the case with (without) 50(5) symmetry. 
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